Controlling Chaotic Transport on Periodic Surfaces 
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We uncover and characterize different chaotic transport scenarios on perfect periodic surfaces by 
controlling the chaotic dynamics of particles subjected to periodic external forces in the absence of 
a ratchet effect. After identifying relevant symmetries of chaotic solutions, analytical estimates in 
parameter space for the occurrence of different transport scenarios are provided and confirmed by 
numerical simulations. These scenarios are highly sensitive to variations of the system's asymmetry 
parameters, including the eccentricity of the periodic surface and the direction of dc and ac forces, 
which could be useful for particle sorting purposes in those cases where chaos is unavoidable. 
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Introduction.— Controlling the transport of particles 
on periodic potential energy surfaces is an old and 
ubiquitous problem appearing in different fields such as 
physics, chemistry, and biology ill]. Specific examples 
include colloidal transport in arrays of optical tweezers 
0, flux creep through type-II superconductors and 
Bose- Einstein condensates with periodic pinning sites [1] , 
among many other. Previous theoretical analysis of the 
motion of particles on surfaces joT-floj considered meso- 
scopic models owing to the great complexity of the differ- 
ent transport scenarios. While non-chaotic regimes have 
been wi dely studied in the context of noisy overdamped 
models [11| and the chaotic regime has been mainly con- 
sidered when directed transport is induced by symmetry 
breaking [t| [l2j], to the best of our knowledge, the fun- 
damental case of deterministic chaotic transport in the 
absence of a ratchet effect has not been considered in 
detail as yet. The study of such a chaotic transport on 
simple periodic surfaces could indeed shed some light on 
diverse chaotic phenomena of great complexity appear- 
ing for example in magnetotransport on antidot lattices 



Model— In this Letter, we consider the classical dy- 
namics of a dissipative particle moving on a standard 
separable periodic potential, with an external force hav- 
ing both dc and ac components, and neglecting thermal 
effects: mx + dV/dx — —fix + f cos9 + f\ x cos (uj x t), 
my + dV/dy — —fiy + fo sin 9 + f\ y cos (u) y t), where 
an overdot denotes a derivative with respect to t, 9 
describes the direction of the dc force fo, /i is the 
phenomenological coefficient of friction, and V(x, y) = 
Vq [cos (2ttx/X x ) + cos (2iry/X y )} /2 is the potential with 
X Xl X y being the characteristic length scales. A main pur- 
pose of the present work is a theoretical characterization 
of the different chaotic transport (CT) scenarios by pro- 
viding analytical estimates of the threshold conditions in 
parameter space by using Melnikov analysis (MA). For 
the sake of a dimensionless description, we put the equa- 



tions of motion into the form 



r x + sin 7-3 = -jr x 
sin r„ 



F 0x cos 9 + F lx cos (Q x t) , (1) 
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where all variables and parameters are dimensionless, 
an overdot denotes a derivative with respect to r = 

1/2 

7rf(2Vb/m) /X x , r x = 2irx/X x ±ir, r y = 2iry/X y ±ir, 
7 = fiX x {2mV ) /it, F 0x = f X x / (nV ), F lx = 

1 /2 

X x fi x / (Wo), Q x = u x X x [m/ (2Vb)] /it, a = X y /X x , 



fly/ fix, and 



Wy/ Ul x 



It is also assumed 



that the system [Eqs. (l)-(2)] satisfies the MA require- 
ments, i.e., the dissipation and forcing terms are small- 
amplitude perturbations of the underlying conservative 



pendula r x>y + sinr^ = (see [l4 - |l6l | for general back- 
ground). Straightforward application of MA to Eqs. (1) 
and (2) yields the Melnikov functions (MFs) 



M£ (to) 



Dt ± 27rFizSech ( ^ ) cos (fi x r ) , (3) 



M±(r ) = ±2TrabF lx sech 



{ 7raci7 x 



cos (cfi x rojj 



respectively, where the positive (negative) sign refers 
to the top (bottom) homoclinic orbit of the conserva- 
tive pendulum, and = ±2ttFq x cos 9 — 87, = 
±2waFo x sin 9 — 807. Since the MFs (3) and (4) have an 
infinity of simple zeros, a main conclusion is that neces- 
sary conditions for the onset of chaotic instabilities are, 
respectively, 



Flx > gfriigyg-iW"* 



2tt 



V 2 



(•5) 



min{|£>+| , \D |} / iracn x \ 
Flx> 2^6 C ° S H^J- (6) 

Next, one can compare the theoretical predictions and 
Lyapunov exponent (LE) calculations [15| with the 
caveat that one cannot expect too good a quantita- 
tive agreement between the two kinds of results because 
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LE provides information concerning solely steady chaos, 
while MM is a perturbative method generally related to 
transient chaos [Tf3]. To quantify the sorting capabil- 
ity associated with the threshold of chaotic transport, 
we evaluate the Cartesian components of the velocity, 
(vi) = lim-r-^oo (fiir)) Jt (i = x,y), where brackets in- 
dicate average over initial conditions, and construct the 
velocity components parallel and perpendicular to the 
external dc force fo, (vn) — (v x ) cos9 + a (v y ) sin#, and 
= — ( v x) sin9 + a (v y ) cos0, respectively. We char- 
acterize the deviation of (v) from f by means of the 
quantifier 

tana = (7) 

where a is the deflection angle [Icj. For the sake of 
clarity, we shall consider here the case with equal fre- 
quencies (c = 1) and both dc and ac forces acting in 
the same direction (f\ x = f\COs9, f\ y = /isin0, and 
hence b — tan 9). By defining Fx = X x fi/ (ttVq), one has 
F\ x = Fx cosO and hence Eqs. (5) and (6) reduce to 



Ft > Ff 



th 



^n{\Dt\,\D-\} 

2tt|cos6»| V 2 



l ' th 2m sin0 



7T0 



n, 



(8) 
(9) 



respectively, where Ff th ,Ff th are the chaotic threshold 
amplitudes. 

Symmetry analysis.— Equations (8) and (9) tell us that 
the onset of chaos in both directions strongly depends 
upon the external force direction 9, which can thus be 
used as a /iig/i-sensitivity control parameter to suppress 
and strength CT in one or the another direction at will. 
Specifically, one straightforwardly obtains from Eqs. (8) 
and (9) that the chaotic threshold amplitudes exhibit (as 
functions of 9) the symmetries: 



F? tth (tt/2±0) 
Fl th (n/2±0) 
Ff >th (tt/4±0) 
F? th (3tt/4±0) 



Ff tth (n/2 T 8) , (10) 

Kth (V2 T 0) , (11) 

KthW*T6), (12) 

Fl th (3n/4 T 0). (13) 

Now, the following remarks are in order. First, Eqs. (10) 
and (11) are valid for any spatial potential (a > 0), while 
Eqs. (12) and (13) are solely valid for a symmetric po- 
tential (a = 1). Second, symmetries (12) and (13) imply 
that different transport regimes are expected in the x- 
and ^-directions as the external force direction deviates 
from the "symmetric" angles 7r/4 and 37r/4, respectively. 
Third, in the absence of multistability (i.e., when a sin- 
gle attractor exists for all initial conditions), symmetries 
(12) and (13) also imply (v x ) (tt/4 ± 6) ~ (v y ) (tt/4 =f 0) 
and (v x ) (37r/4 ± 9) ~ (v y ) (37r/4 =p 9), respectively, and 
hence tana (as a function of 9) exhibits the symmetry 



tana (7r/4 - 
tana (37r/4 - 



tana (tt/4 — 
tana (37r/4 



(14) 
(15) 



i.e., for a symmetric potential, tana is an odd function of 

8 with respect to the angles 7r/4 and 37r/4, respectively. 
Note that this is no longer the case for an asymmetric 
potential according to the first remark. 

Numerical results.— Extensive numerical simulations 
confirmed all the above theoretical predictions. Thus, 
by varying 9 one can find different transport regimes 
(see Fig. 1, top panel): CT in both directions (as 
for 9 = {27r/9,5vr/18}), CT in one direction while 
intermittent periodic transport (PT) in the other 
(as for 9 = {71-/6, 7r/3}), PT in both directions (as 
for 9 = {77r/36,7r/4,ll7r/36}), and PT in one di- 
rection while periodic oscillation in the other (as for 

9 = {137r/36, 57r/36}). Since the onset of chaos also 
depends upon the particle mass (through the coefficient 
of friction, cf. Eqs. (8) and (9)), such an 6*-dependcncc 
can therefore be used to sort different particles according 
to their mass. For two kinds of particles with different 
masses, this means that one can obtain analytical 
estimates of the optimal force directions, 9 opt , from Eqs. 
(8) and (9) such that one particle exhibits CT while the 
other does not, the remaining parameters being held con- 
stant. Numerical experiments confirmed this scenario as 
is shown in Figs. 1 (medium panel) and 2. Additionally, 
the onset of chaos also depends upon the eccentricity 
parameter a (Eq. (9)): Decreasing or increasing a from 
1 (symmetric potential) means increasing the potential's 
asymmetry. Thus, the eccentricity of the periodic 
potential can also be used as an effective parameter to 
control CT on a periodic surface, as in the case of optical 
potentials for example [l7j . Figure 1 (bottom panel) 
shows an illustrative example where typical trajectories 
are plotted for increasing values of a from 1. Starting at 
a situation where CT occurs in both directions (a = 1), 
one finds that increasing the potential's asymmetry 
(a > 1) changes the motion to PT in both directions (as 
for a = 1.2). This behaviour changes again to CT in 
both directions for higher values of a (as for a = 1.4), 
and finally changes to PT in the x-direction while remain 
bounded inside a well in the y-direction (as for a = 1.6). 
Also, the effectiveness of a fixed external force at sorting 
heavy particles is enhanced by breaking the potential 
symmetry (recall that 7 ~ mT 1 / 2 , see Fig. 2). Figure 
3 shows illustrative instances of maximal LEs, A+ and 
A+, which quantify the chaotic dynamics in the x- and 
y-directions, respectively, versus 9 for two values of the 
eccentricity parameter. Remarkably, these diagrams 
present relevant symmetries which are coherent with 
those of the chaotic threshold amplitudes [Eqs. (10)- 
(13), respectively]: A+ (tt/2 ± 9) = A+ (n/2 =p 0), 
A+(n/2±9) = A+ (n/2 T 0), A+ (tt/4 ± 9) = 
A+(tt/4t0), A+(3tt/4±60 = A+ (3tt/4 t 0). It is 
worth mentioning that this coherence is far from trivial 
in the sense that, to the best of our knowledge, there is no 
theoretical connection between MA predictions and LEs 
for the present system, thus indicating the relevance and 
depth of the chaotic threshold symmetries in parameter 
space. One typically finds how different chaotic and non- 
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FIG. 2: Deflection angle vs coefficient of friction for two 
values of the eccentricity parameter: a = 1 (0) , 1.5 (A). The 
inset shows the deflection angle vs eccentricity parameter for 
7 = 0.2. Other parameters are: 9 = tv/5,Fq x = 0.28, Fi = 
1,Q X = 0.68. The solid lines are solely plotted to guide the 
eye. 




FIG. 3: Maximal LEs A+, A+ as a function of the angle 9 for 
two values of the eccentricity parameter: a — 1 (top panel) , 
1.5 (bottom panel). Other parameters are: Fox = 0.28, Fi = 
1,7 = 0.25 and Q x = 0.68. 



FIG. 1: Top panel: Trajectories for a net force applied at 
different angles: 9 = 13tt/36 (a), vr/3 (b), llvr/36 (c), 5vr/18 
(d), tt/4 (e), 2tt/9 (f), 7tt/36 (g), vr/6 (h), and 5tt/36 (i) for 
a — 1 and 7 = 0.25. Medium panel: Trajectories for a — 
1,9 = 27r/9, and four values of the dimensionless coefficient 
of friction: 7 = 0.25 (a), 0.3 (b), 0.35 (c), and 0.4 (d). Bottom 
panel: Trajectories for 9 — 27r/9,7 = 0.25, and four values 
of the eccentricity parameter: a — 1 (a), 1.2 (b), 1.4 (c), 
and 1.6 (d). Other parameters are: Fo x = 0.28, Fi = 1, and 
fix = 0.68. Dotted lines indicate the direction of the external 
force. 



chaotic regimes drastically change over certain 9 ranges 
as the potential becomes asymmetric. For instance, PT 
in both directions at 9 = {7r/4,37r/4} for a symmetric 
potential (a = 1) changes to CT in solely one direction 
for an asymmetric potential (a = 1.5) (cf. Fig. 3). 
Finally, numerical simulations confirmed the accuracy of 
predictions (14) and (15) as is shown in Fig. 4. Starting 
with CT in both directions at 9 = 7r/4 for a symmetric 
potential (Fig. 4, top panel), one sees that the deflection 
of particles increases as 9 deviates from ir/4 according 
to the route described in Fig. 1, top panel. Maximum 
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deflection occurs at symmetric angles with respect to 
t/4, ^Tx (t/4 - - C ~ t/4) , where there 
is PT in one direction while periodic oscillation in 
the other. For 6 ^ 6>™| x (0 sC 6»^ x ) , this transport 
regime remains, i.e., (v x ) = {{v y ) = 0) and hence 
tan a (OCT*) = cot 6 (tan a (6 < fe) = - tan 6) 
(cf. Eq. (7)). For an asymmetric potential (Fig. 4, 
bottom panel), the dependence of the deflection angle on 
the external force direction essentially presents a similar 
scenario to that of the symmetric case, but now tan a is 
no longer an odd function with respect 7r/4, as predicted 
(cf. third remark). 
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FIG. 4: Deflection angle vs external force direction for 
Fox = 0.28, Fi = 1, fl x = 0.68, two values of the eccentricity 
parameter: a — 1 (top panel), 1.5 (bottom panel), and differ- 
ent values of the coefficient of friction: 7 = 0.25 (0) , 0.4 (a) 
(top panel) and 7 = 0.15 (0) , 0.25 (A) (bottom panel). Also 
plotted are the functions — tan 6 and cot 9 (dashed lines, see 
the text ). 



Conclusions.— To sum, we have demonstrated theoret- 
ically and numerically through a simple and general sys- 
tem that reliable control of sorting on periodic surfaces 
is achieved for chaotic particles by identifying the rele- 
vant symmetries of the chaotic threshold in parameter 
space. We uncovered and characterized different sort- 
ing scenarios associated with symmetric and asymmetric 
spatial potentials, which could motivate experiments in 
different contexts such as optical and antidot lattices. 
Among the most interesting extensions of this work are 
the case with the ac and dc forces having different direc- 
tions, where preliminary results indicate the presence of 
intriguing "absolute negative mobility" phenomena [l8j . 
as well as the study of the effect of noise on the present 
transport scenarios: Even very small amounts of noise 
may cause both a transition from a bounded state to 
a running state and a significative modification of the 
chaotic threshold in parameter space [l9| . Our current 
work is aimed at exploring these cases. 
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